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1. Introduction 

The inclusion of the Coulomb interaction in the description of the three-particle 
scattering is a challenging task in theoretical few-body nuclear physics. Due to its 
long range, the Coulomb potential w{r) does not satisfy the mathematical prop- 
erties required for the formulation of the standard scattering theory. There is a 
number of suggestions how to ove rcome t his difficulty; most of them are based on 
the configuration-space frameworl(P2l3ll] and are limite d to energies below three- 
body breakup threshold (3BBT), while the otherJSEEl have not matured yet into 
practical applications. Up to now only few approaches led to the results above 
3BBT. Those are configuration-space calculations for proton-deuteron [p-d) elastic 
scattering using the Kohn variational principlJ^l and the screening and re nqrrnaliza- 
tion approach in the framework of momentum-space integral equations j^^^^^^^^^^ 
the latter method will be discussed in more details. Very recently p-d results above 
3BBT were also obtained using modified Faddeev equation in configuration space 
together with the dumping (screening) of particular Coulomb contributions.^^ 



2. Method of screening and renormalization 

Our treatment of the Coulomb interaction is based on the idea of screening and 
renormalization proposed in Refs.[T31[TS]for the scattering of two charged particles. 
The standard scattering theory is formally applicable to the screened Coulomb 

1 



October 24, 2008 15:22 WSPC/INSTRUCTION FILE deltuva 



2 A. Deltuva 

potential which, in the r-space representation, we choose as 

WR{r) ^w{r) e-^''/^^", (1) 

where R is the screening radius, and n controls the smoothness of the screen- 
ins. In 1974 TayloJl^l suggested that even for the description of systems with the 
Coulomb interaction the standard scattering theory may be useful, since in nature 
the Coulomb potential is always screened. The study of the two-particle system 
with the screened Coulomb interaction revealed that, as expected, the physical ob- 
servables become insensitive to screening provided it takes place at sufficiently large 
distances R and, in the i? — *■ oo limit, coincide with the corresponding quantities 
known from the analytical solution of the two-particle Coulomb problem: though 
the on-shell screened Coulomb transition matrix (p/|t_R(ei -|-iO)|pi) with energy 
and momenta Pf = Pi diverges in the i? — > oo limit, after renormalization by (an 
equally) diverging phase factor z^^(pi) it converges as a distribution to the well 
known proper Coulomb amplitude {pf\tc\pi), i.e., 

lim {pf\tR{e^ + iO)\p^)z]^'^{p,) = {pf\tc\Pt). (2) 

R — 'oo 
_l 

Renormalization by Zj^^{pi) in the R ^ oo limit relates also the screened and 
proper Coulomb wave functions.^l^ The convergence of the scattering amplitude as 
a distribution corresponds to the physical conditions in a real experiment where the 
incoming beam is not a plane wave but a wave packet and the forward scattering 
cannot be measured.!^ This justifies the replacement ([2]) in practical calculations. 
We emphasize that the renormalization ([2|) relates only the scattering amplitudes; 
its application to the off-shell transition matriceJl^ in unjustified. However, for the 
calculation of observables, only the renormalization of the on-shell transition matrix 
and wave function is needed. 

The scree ning and renormalization approach can be applied to more compli- 
cated systems.^^ Here we briefly recall the procedure which is described in details 
in Refs.[TTl[T2l In the multichannel on-shell transition matrix {f\Tlj^\Ei + iO)\i) be- 
tween initial and final channel states \i) and |/), Ef — Ei, derived from nuclear plus 
screened Coulomb potentials, one has to isolate the diverging screened Coulomb con- 
tributions in the form of a two-body on-shell transition matrix and two-body wave 
function with known renormalization properties. This can be achieved using the two- 
potential formalism as long as in the initial/final states there are no more than two 
charged bodies (clusters). This also enables to decompose {f\Tp^ {Ei -I- into 
contributions with different range properties. The long-range part, the two-body 
on-shell transition matrix {f\T^ji"{Ei -\- «0)|j), derived from the screened Coulomb 
potential of the form ([T]) between the centers of mass (cm.) of the two charged 
bodies in the initial state, is present in the elastic scattering only. After renormal- 
ization by the diverging phase factor Z~j^, this contribution converges towards its 
i? — > cxD limit very slowly, and, in general, as a distribution only, but the result 
(/l^ac' N)' pure Coulomb amplitude of two-body nature, is known analytically. 
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The remaining part of the elastic scattering ampHtude as well as the amplitudes 
for transfer and breakup are short-range operators that are externally distorted by 
Coulomb. Due to their short-range nature, convergence with R after the renormal- 
ization by the corresponding phase factors is fast and, therefore, the i? — > oo limit 

can be calculated numerically with high accuracy at finite R. Thus, the physical 

(a) 

scattering amplitudes are obtained after rcnormalization of {f\TpJ{Ei + iO)|i) in 
the i? oo limit as 

+ hm {Z]i (/I [T^f ) [E, + zO) - 5p^T-^il- (E, + zO)] |z)Zr/ }. 

One can use standard scattering theory to calculate the multichannel transition 
operators {Ei -\- iQ) at finite screening radius R and make sure that R is large 
enough for the convergence of the results. We solve Alt, Grassberger, and Sand- 
has (AGS) equations for three- and four-particle scatterinim 

which are equiv- 
alent to Faddeev and Yakubovsky equations .^^ll^ We employ momentum-space 
partial- wave representation as described in detail in Refs. [23l [24l |25]for three- and 
four-nucleon scattering without the Coulomb force. However, the screened Coulomb 
interaction, due to its longer range, compared to the nuclear interaction, brings ad- 
ditional difficulties: quasisingular nature of the potential and slow convergence of 
the partial-wave expansion. The right choice of the parameter n controlling the 
smoothness of the screening ([T]) is essential in resolving those difficulties as we 
demonstrated in Rcf. 1111 We work with a sharper screening than the Yukawa screen- 
ing (n = 1) of Refs. iQl ITOl We want to ensure that the screened Coulomb potential 
wr{t) approximates well the true Coulomb one w{r) for distances r < R and simul- 
taneously vanishes rapidly for r > R, providing a comparatively fast convergence of 
the partial- wave expansion and less pronounced quasisingularities. In contrast, the 
sharp cutoff {n oo) yields an unpleasant oscillatory behavior in the momentum- 
space representation, leading to convergence problems. Depending on the reaction 
we found the values 3 < n < 8 to provide a sufficiently smooth, but at the same 
time a sufficiently rapid screening around r = i?. In any case the screening radius R 
needed for convergence in Eq. ^ is considerably larger than the range of the nuclear 
interaction and, therefore, the calculation of TpJ{Ei -\- iQ) requires the inclusion of 
partial waves with angular momentum much higher than needed for the nuclear 
potential alone. This problem can be solved in an efScient and reliable way either 
by using the perturbative approach for high two-particle partial waves, developed 
in Ref. |26l or even without it as discussed in Ref. '27', 

The internal criterion for the reliability of our method is the convergence of 
the observables with screening radius R used to calculate the Coulomb-distorted 
short-range part of the amplitudes in Eq. ^ . Numerous examples for three-nucleon 
hadronic and electromagnetic reactions, a-d and p-'^He scattering can be found in 
Refs. [TTl[T^[^[^ [5D. In most cases the convergence is impressively fast; the screen- 
ing radius i? = 10 to 30 fm is sufficient. The exceptions requiring larger screening 
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Fig. 1. ^So phase shifts of pp scattering at 0.01, 0.1 and 1 MeV cm. energy. Results obtained 
using screening and renormalization method are shown as function of the screening radius R (solid 
curves). Dashed lines are the exact results. 



radii are the observables at very low energies and the breakup differential cross sec- 
tion in kincmatical situations characterized by very low relative energy Ej-d between 
the two charged particles, e.g., p-d breakup or photodisintegration of ■^He close to 
the pp final-state interaction (pp-FSl) regime.'^ In there, the Coulomb repulsion 
is responsible for decreasing the cross section, converting the FSI peak obtained in 
the absence of Coulomb into a minimum with ze ro cros s section at Ei-c\ ~ 0. Such a 
behavior is seen in the experimental data as well .1^1122] ^pj-^g glow convergence under 
those conditions is not surprising, since the renormalization factor itself as well as 
the Coulomb parameter become ill-defined, indicating that the screening and renor- 
malization procedure cannot be applied at E^d = 0. Therefore an extrapolation has 
to be used to calculate the observables at E'rci — 0, which works pretty well since 
the observables vary s moo thly with i?rci- 

It has been claimecJI^ that the method of screening and renormalization ceases 
to work even for pp scattering below 0.1 MeV cm. energy. In Fig. [1] we prove that 
this is not so by getting well converged results at 0.01 MeV. However, Fig. [T]also 
shows that the screening radius indeed has to be increased with decreased energy 
and at some point would become too large for reliable numerical calculation. 



3. Results and summary 

We have discussed how the Coulomb interaction between the charged particles can 
be included into the description of few-body reactions using the old idea of screening 
and reno rmahz ation.ll^ The calculations are done in the framework of AGS integral 
equationJISI^S! in momentum-space. T he scree ning and renormalization approach 
has already been used for p-d scatterinJ^HnESl but with limited success: those cal- 
culations were based on quasiparticle equations with rank-1 separable potentials 
and, in addition, the screened Coulomb transition matrix was approximated by the 
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screened Coulomb potential; none of these approximations is used by us. It is the 
new screening function that allows us to avoid these approximations and obtain fully 
converged results. Furthermore, the results for p-d elastic scattering obtained by the 
present technique were comparetJ^^l with those of Ref. [8] obtained from the varia- 
tional solution of the three-nucleon Schrodinger equation in configuration space with 
the inclusion of an unscreened Coulomb potential between the protons and imposing 
the proper Coulomb boundary conditions explicitly. Good agreement over a wide 
energy range was found indicating that both techniques for including the Coulomb 
interaction are reliable. At very low energies the coordinate-space treatments remain 
favored since there the method of screening and renormalization converges slowly 
and therefore becomes technically too demanding, but at higher energies and for 
three-body breakup reactions it is more efficient. A detailed comparison still has to 
be done but it seems that there is a reasonable agreeme nt b etween momentum- and 
coordinate-space results also in the case of p-d breakujJISl and p-^He scattering.!^ 
The present method was used to study p-d elastic scattering and breakup in 
Refs. [TTJ [T2J |M1 IMl p-d radiative capture in Refs. [TTJ [37| and photo- and elec- 
trodisintegration of '^He in Refs. [TTl [T^ [351 Furthermore, it was applied to the 
nuclear reactions dominated by three-body degrees of freedom, i.e., deuteron scat- 
tering on stable nuclei or proton scattering on one -ne utron halo nuclei. Examples 
are low-energy a-d elastic scattering and_breakujJ^ and d -\- ^^C and p -\- ^^Be 
elastic, transfer, and breakup reactioni^SEQllT] .^^here also the accuracy of tradi- 
tional approximate nuclear reaction approaches like Continuum Discretized Coupled 
Channels (CDCC) method, Glauber, and Distorted Wave Impulse Approximation 
(DWIA) could be tested. Finally, all elastic and transfer four-nucleon reactions be- 
low three-body breakup threshold, i.e., n -\- "^H, p -\- ■^He, p + ^H, n + ■^He, and 
d-\-d, have been studied in Refs. [25l [30l [42l |43l The conclusion is that the Coulomb 
effect is important at low energies for all kinematic regimes, but gets confined to 
the forward direction in elastic scattering at higher energies. In p-d breakup and in 
three-body e.m. disintegration of ^He the Coulomb effect is extremely important 
in kinematical regimes close to pp-FSI. There the pp repulsion converts the pp- FSI 
peak obtained in the absence of Coulomb into a minimum with zero cross section.!^ 
This significant change of the cross section behavior has important consequences in 
nearby configuration s w here one may observe instead an increase of the cross sec- 
tion due to Coulomb.l^ However, some of the long-standing discrepancies between 
experiment and theory like the space star a nomaly in p-d breakup are not resolved 
by the inclusion of the Coulomb interaction.^S^IBl ^ very strong Coulomb effect is 
found in a-d breakup where the shift of ap P-wave resonance position leads to the 
corresponding shifts of the differential cross section peaks.-^ 
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